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Maintaining quantum coherence is a crucial requirement for quantum computation; hence protect- 
ing quantum systems against their irreversible corruption due to environmental noise is an important 
open problem. Dynamical decoupling (DD) is an effective method for reducing decoherence with a 
low control overhead. It also plays an important role in quantum metrology, where for instance it 
is employed in multiparameter estimation. While a sequence of equidistant control pulses (CPMG) 
has been ubiquitously used for decoupling, Uhrig recently proposed that a non-equidistant pulse 
sequence (UDD) may enhance DD performance, especially for systems where the spectral density 
of the environment has a sharp frequency cutoff. On the other hand, equidistant sequences out- 
perform UDD for soft cutoffs. The relative advantage provided by UDD for intermediate regimes 
is not clear. In this paper, we analyze the relative DD performance in this regime experimentally, 
using solid-state nuclear magnetic resonance. Our system-qubits are 13 C nuclear spins and the en- 
vironment consists of a H nuclear spin-bath whose spectral density is close to a normal (Gaussian) 
distribution. We find that in the presence of such a bath, the CPMG sequence outperforms the UDD 
sequence. An analogy between dynamical decoupling and interference effects in optics provides an 
intuitive explanation as to why the CPMG sequence performs superior to any non-equidistant DD 
sequence in the presence of this kind of environmental noise. 

PACS numbers: 03.65.Yz,03.67.Pp,76.60.-k ,76.60.Lz 



I. INTRODUCTION 



Quantum information processing (QIP) relies on stor- 
ing and manipulating information in quantum mechan- 
ical states associated with accessible subunits called 
qubits For reliable QIP, the information encoded in 
the quantum register (the collection of qubits) must be 
retained for arbitrarily long time [2]. However, in any 
physically realizable QIP architecture, the qubits can- 
not be completely isolated, but are weakly coupled to 
a large number of degrees of freedom of their environ- 
ment (bath) . This causes the corruption of the quantum 
state associated with the qubit, a process known as deco- 
herence H . This process limits the time scale over which 
quantum information can be retained j| and the distance 
over which it can be transmitted 0"l2|- 

Combating this decoherence process to extend the life 
time of quantum states or processes and the distance 
bounds to transmit them is a necessary step in building 
a quantum computer One of the simplest and most 
effective techniques suggested for this purpose is Dynam- 
ical Decoupling (DD) [SEl- 

It also has several applica- 
tions in quantum metrology, for example in the magne- 
tometry using single spins |lll4l3j and multiparameter es- 
timation [T3|. DD consists of the application of 7r-pulses 
to the qubits, which revert the decay due to the system- 
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environment interaction. The simplest implementation 
of this method is the Hahn echo experiment [lj| . It is 
employed in the case when the system-environment (SE) 
interaction is a pure dephasing process - i.e. there is no 
net exchange of energy between the system and the bath. 

However, it is immediately obvious that any compo- 
nent of the interaction with the bath that varies on a 
time scale shorter than r, the delay between the pulses, 
cannot be refocused [HI, [l?} . Hence the effectiveness of 
the Hahn echo and its generalization as a train of equidis- 
tant pulses - the CP [HI and CPMG [I?} sequences - de- 
pend crucially on keeping the delays r between the pulses 
sufficiently short. However, this delay is always limited 
by hardware or the maximum power deposition in the 
sample. Moreover, rapid DD pulsing interferes with the 
controls that are necessary for computation [IH, [lj| . For 
example, if one wants to choreograph the dynamics of 
the spins in a particular manner while still requiring that 
coherence survives for long, the requisite controls have to 
be applied in parallel with the DD pulses. Although the 
recent proposal of Dynamically Corrected Gates (DCGs) 
[20M23I ] seeks to address this problem, when r is very 
short there is an enormous demand on the hardware. 

Hence there is a strong motivation for finding a DD 
protocol that, for a given number of pulses, provides the 
best performance. For the case of a purely dephasing SE 
interaction, a possible approach to this problem consists 
in visualizing the DD pulses as genera ting a f ilter for the 
SE interaction (environmental modes) [24j|25j . For exam- 
ple, the CPMG sequence acts as a band-stop or high-pass 
filter (with bandwidth ui = 2it/t) since any interaction 
component that varies slower than r is refocused and hi- 
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tered out. The resulting decay rate of the quantum state 
is determined by the overlap of the spectral density of 
the bath with the filter generated by the pulse sequence 
[24l [25j . If one knows the exact form of this function, 
it may be possible to design a suitable filter, and hence 
a DD sequence, that leads to the slowest decay {25U31I ]. 
This was the motivation behind the DD sequence sug- 
gested by Uhrig (2(| . It provides the non-equidistant dis- 
tribution of pulses that causes the flattest band-stop filter 
around w = [H|,[32| , and is considered the best sequence 
for combating low frequency noise. Experimental work 
[HiUlil showed that indeed UDD outperforms CPMG 
for certain types of noise with a sharp high-frequency cut 
off. 

However, when the spectral density of the bath has a 
long tail (soft cutoff), the CPMG sequence has been pre- 
dicted [II HI and found to outperform the UDD 
sequence. In an intermediate regime, for example a spin- 
bath where the spectral density is Gaussian [371 ] . it was 
shown rece ntly that CPMG outperforms the UDD se- 
quence [38H41| . However, a rigorous comparison of both 
these sequences in such regimes is a matter of current re- 
search [H, HIHm . Of particular interest is determining 
the conditions (for example, the maximum allowable r) 
under which the UDD continues to provide an advantage 
over CPMG. 

This was tackled in a purely dephasing bosonic bath by 
Hodgson et at [43J. They predicted that DD protocols 



with non-equidistant pulses like UDD lose their advan- 
tages when lower bounds exist for the pulse separation r. 
Uhrig and Lidar j46j presented analytical performance 
bounds for UDD. They considered instantaneous per- 
fect pulses and bounded environment and generic system- 
environment Hamiltonians. They showed that for a fixed 
total duration, the survival probability can be increased 
by increasing the number of pulses by incrementing the 
DD order. However, if a minimum delay between pulses is 
imposed, because the sequence time scales with the num- 
ber of pulses (DD order) , they predict that longer cycles 
(higher orders) are not always advantageous. These pre- 
dictions are based on perturbative treatments of the SE 
interaction. Accordingly, they apply to the fast control 
regime where the cycle time is sufficiently short for the 
perturbative treatment. Because these approaches can- 
not lead to a lower bound on the attainable DD error in 
the presence of timing constraints, recently Khodjastch 
et al. UK obtained bounds for the UDD performance as 
a function of its order by using an alternative method for 
bosonic baths. 

In light of this, it is important to provide a fair and 
practically relevant comparison between these different 
DD approaches. Here, we use the benchmarking method- 
ology of our previous work [38| : which DD sequence max- 
imizes the survival time of the quantum coherence for a 
given average spacing between the pulses? In the earlier 
paper, we concentrated on a comparison between differ- 
ent sequences with equidistant pulses; here, we analyze 
the effect of a variable pulse spacing, using the same 



approach. The spectral density of our bath is roughly 
Gaussian, which may be considered to be intermediate 
between a sharp cutoff and a long tail. 

Several recent experimental results [381441} indicate 
that CPMG performs better than UDD when the qubit 
is coupled to a spin-bath. The superior CPMG perfor- 
mance has been attributed to the fact that it compen- 
sates pulse errors along the direction of the rotation axis 
41|. while such compensation is reduced in UDD 
331 . Also it was attributed to a possible soft-cut 
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off of their model system [39l - l4l| . Our experiments and 
simulations indicate that even in the case of ideal pulses, 
the CPMG outperforms UDD. In fact, we conjecture that 
for a purely dephasing SE interaction, when the spectral 
density is a Gaussian distribution and a rapidly fluctu- 
ating environment, the CPMG sequence would maintain 
the state of the qubit system for the longest period. Our 
reasoning is based on determining how the DD pulse se- 
quence alters the time-averaged SE interaction that the 
system "sees", for each frequency component of this inter- 
action [24|. Drawing an analogy with optics, each pulse 
of the Uhrig sequence modifies the SE interaction such 
that the interactions in different windows "interfere" de- 
structively with each other. This destructive interference 
is perfect for slow interactions. While the UDD protocol 
provides the flattest stop- band at uj = (24| , the CPMG 
protocol provides the widest. When a DD sequence is 
applied repetitively (say M times) , it becomes analogous 
to passing the different frequency components of the SE 
interaction through a diffraction grating. Overall, the 
optical analogy presented in this paper provides an in- 
tuitive way to understand how the shape of the spectral 
density function affects DD performance. 

This paper is organized as follows. Sec. HJdcscribes the 
physical system that we use in our experiments. In Sec. 
IIIII we review the key ideas of dynamical decoupling. Sec. 
IIVI provides the pulse sequence protocols of UDD |26l.l32l| 
and CPMG QJ, [13 ■ Sec. [V] contains the results of our 
experiments. Sec. IVII describes a semi-classical model 
that explains the experimental results, and suggests that 
the CPMG sequence is the best DD sequence for com- 
bating a purely-dephasing SE interaction in a spin-bath. 
Within this section we draw an optical analogy to DD 
sequences. 



II. THE SYSTEM 

For our experimental evaluation of the relative perfor- 
mance of DD sequences, we use the simplest possible sys- 
tem consisting of a single qubit and a purely dephasing 
environment consisting of a spin-bath. For the system 
qubit, we use 13 C nuclear spins (S = 1/2), and for the 
environment 1 H nuclear spins (/ = 1/2) that act as the 
spin-bath. We use polycrystalline adamantane where the 
natural abundance of the 13 C nuclei is about 1%, and to a 
good approximation each 13 C nuclear spin is surrounded 
by about 100 1 H nuclear spins. The 13 C- 13 C interaction 
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Figure 1: Schematic representation of dynamical decoupling 
The solid boxes represents the control pulses. 



is completely negligible compared to the 13 C- 1 !! and 1 H- 
*H interactions. The spin-spin interaction is dominated 
by the dipolar interaction [37] . 

We shall refer to the spin operator of the system qubit 
( 13 C) as S, and to the j th bath spin ( 1 H) as P. The 
Zeeman frequencies of the system- and bath spins are uis 
and w/, respectively. The dipolar coupling constants for 
the system-bath interaction are bsj, and for the intra- 
bath interaction dij. 

It is convenient to describe the dynamics of the system 
in a rotating frame of reference [371 ] . where the system 
rotates at the (angular) frequency ujs around the z-axis 
and the environment at wj. The total (free evolution) 
Hamiltonian is then 



H f =Hi 



T-Lse + He, 



(1) 



where T-Ls is the system Hamiltonian, He the environ- 
ment Hamiltonian and Use the system-environment in- 
teraction Hamiltonian: 
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This type of system is encountered in a wide range 
of solid-state spin systems, as for example electron spins 
in diamonds [39|, |4l|, [47] , electro n sp ins in quantum dots 
[iol l48l |49| and donors in silicon [Hof, [H| , which appear to 
be promising candidates for future QIP implementations. 
In particular, we consider the case where the interaction 
with the bath is weak or comparable with the intra-bath 
interaction. 



III. DYNAMICAL DECOUPLING 

Dynamical decoupling aims to reduce the interaction of 
the system with the environment by applying sequences 
of short, strong pulses that invert the system qubits 
!). lOj. We write Hc(S){t) f° r tnc corresponding con- 
trol Hamiltonians. It is assumed that the environment 
cannot be directly controlled. 
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Figure 2: (Color online) Distribution of pulses for (a) CPMG 
and (b) UDDio sequences. Here r is the average distance 
between pulses. The arrows denote the length of a single 
cycle of period t c . Note that a CPMG cycle consists of 2 
pulses while an N th order UDD cycle has N pulses. 



DD sequences usually consist of cycles of pulses. Fig- 
ure [1] shows an example of such a cycle. In the rotating 
frame, the operator that describes the evolution of the 
total system from to t c , the duration of the cycle, is 



N 



(5) 
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U(r c ) = U f (r N+1 )l[Ub(r P )Uf(T i ) 

i=l 

where the free evolution operator is 

U f (t) = exp{-iH f t} 
and the control evolution operators are 

U l c (r p ) = Texp j-» j " dt< (U f + H c{s){ 



where T is the Dyson time-ordering operator [52], [53] and 
t p the pulse duration. Let U represent the time at which 
the i th control operation starts. Then, the delay times 
between the control Hamiltonians are Tj = £»— (Aj-1 + t p ) 
for i = 2, ..,n + 1 and n = t\ — to, where to = and 
£jv+i = t c . If the basic cycle is iterated M times (see 
Fig. [T]), the total evolution operator becomes 



(*))}, (7) 



U(t = Mt c 



U{t c ) 



M 



(8) 



The propagator ([5]) for a single cycle of duration r c can 
be written in terms of an effective Hamiltonian: 



U(r c ) 



-m. 



(9) 



Using average Hamiltonian theory [54j H e ff can be ex- 
pressed as a series expansion, using, e.g., the Magnus 
expansion (55|, 



ileff = n {0) + n {2) + ... = Y j h 



(m) 



(10) 
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The lowest order term T-L^ is the average over the period 



H(t)dt. 



(11) 



A DD sequence with ideal pulses makes W^> =We , i.e., 
the interaction Hamiltonian vanishes to zeroth order. In 
addition, the DD sequences are designed such that the 
hig her order terms have reduced norm or vanish [ToL l26l . 

As we discussed in the previous section, the 13 C- 13 C 
dipolar interaction is negligible compared with the 13 C- 
: H interaction. However, when DD is active and removes 
the 13 C- 1 H interaction, the 13 C- 13 C couplings become 
the dominant interaction for causing decoherence. This 
holds strictly for ideal pulses, which do not affect the 13 C- 
13 C interaction. Under these conditions, the decay rate 
is limited by the 13 C- 13 C couplings. However, if pulses 
have finite duration, they can further extend the lifetime 
of the component of the magnetization aligned with the 
effective field (EMI. 



IV. PULSE SEQUENCES 

A. Carr-Purcell (CP) and 
Carr-Purcell-Meiboom-Gill (CPMG) 

The Carr-Purcell (CP) sequence [l6| was first intro- 
duced as a means to suppress inhomogeneity of the Bq 
field A£?q seen at the j th nuclear site when the molecules 
containing the spins are undergoing diffusion. The re- 
sulting time-dependent system-environment interaction 
interferes with the refocusing process of a Hahn echo [l5[ 
if the time dependence happens on a time scale compa- 
rable to or faster than the refocusing time. 

In our solid-state spin system, the time dependence of 
the environment arises from fluctuating dipolar fields due 
to the neighboring spins of the bath [54j . This generates 
a time-dependent inhomogeneity. By suppressing it, the 
CP sequence acts as a method of dynamical decoupling 
since, in effect, it suppresses the SE interaction 

The CP sequence consists of a train of refocusing pulses 
generating spin-echoes. If the refocusing pulses are ideal, 
i.e., they cause a complete n rotation of the nuclear 
spins in an infinitesimal time, and the spin-bath is static 
{[HeMse] = 6), then all terms H (m) of the Magnus 
expansion vanish, and the inhomogeneity is completely 
suppressed. However, this does not happen when the 
spin-bath is not static or the pulses have finite duration 
or contain errors. For ideal pulses but a fluctuating en- 
vironment only the average system environment Hamil- 
tonian H^e, the zeroth order of the Magnus expansion, 
is suppressed [38| . 

If flip-angle errors are also considered, the system- 
environment coupling no longer vanishes, even in lowest 
order, T-Lg E ^ 6, and the resulting propagator increases 



with the flip-angle error. Considering this, Meiboom and 
Gill suggested a modification to the CP sequence fl7| . 
now called the CPMG sequence, where the rotation axis 
of the pulses is the same as the orientation of the initial 
state of the spins. In this case, flip-angle errors have no 
effect in zeroth order. Fig. [2^ shows a CPMG sequence. 
Note that the CPMG cycle requires only 2 pulses; hence 
a sequence of 10 pulses covers 5 cycles. In the figure, r 
is the average distance between pulses, while r c denotes 
the length of a single cycle. 



B. Uhrig dynamical decoupling (UDD) 

In the CPMG sequence, the separation between ad- 
jacent pulses is constant throughout the sequence. In 
a seminal paper (26l . |32| . Uhrig relaxed this condition, 
searching for the optimal combination of delays that 
would minimize the effect of a purely dephasing system- 
environment coupling. He showed that the optimal dis- 
tribution for reducing low-frequency noise corresponds to 
a sinusoidal modulation of the pulse delays. More specif- 
ically, one cycle of an N th order UDD sequence consists 
of N pulses applied at times 



U 



t c sin 



2(N+1) 



(12) 



where t^ + i = t c is the cycle time and t = the starting 
time. Fig. [2Jd shows the UDD sequence for N = 10. The 
CPMG sequence is the simplest UDD sequence with or- 
der N — 2. The UDD sequence was shown to be universal 
for any purely dephasing Hamiltonian [fril ] . 



V. EXPERIMENTS 

A. Measurement scheme 

The general procedure used in our experiments is il- 
lustrated in Fig. |3l Cross-polarization from the abun- 
dant proton spins to the 13 C spins is used at the initial 
preparation stage to increase the 13 C polarization [o5| . 
After this transfer, we store the enhanced 13 C polariza- 
tion along the z axis and wait for a time longer than 
the dephasing time to erase any correlation that could 
arise during the cross-polarization process. The magne- 
tization is then rotated to an initial state, represented by 
the density operator po = Si XtV y, in the .Ty-plane that is 
transverse to the strong static field. We use two distinct 
initial states considering the relative phases of the polar- 
ization of po and the pulses to be used for DD (28|, l38j|- 
the "longitudinal" state in Fig. @|a) where the pulses are 
applied in the direction of the po polarization, and the 
"transverse" state in Fig. EJb) where they are applied 
perpendicular to the po polarization. This translates, in 
Fig. El to applying the DD tt pulses with an X phase 
where a tt/2 pulse is applied before of the DD sequence 
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Figure 3: (Color online) The experimental scheme: the initial 
state preparation uses cross-polarization (CP) to enhance the 
polarization of the 13 C nuclear spins. The actual dynamical 
decoupling consists of M cycles of the iV-order DD sequence. 
The DD pulses are applied along the X axis, while the initial 
states are created using a (^) pulse with Y (-X) phase for the 
longitudinal (transverse) initial conditions. 



with a Y phase for the longitudinal initial condition, and 
with a —X phase for the transverse condition. M cy- 
cles of the TV-pulse DD sequence are applied; the dura- 
tion for the DD tt pulses was 10.4 ps. The errors of the 
pulse delays (compared to the theoretical values) are 2.5 
ns on average, limited by the timing resolution of our 
pulse generator. The error of the pulse durations (jit- 
ter) is < 0.1 ns. Therefore timing errors are negligible 
compared to amplitude errors, which are of the order of 
10%. In general for the range of delays between pulses 
used in our experiments, we did not observe substantial 
differences between considering the delays from the pulse 
edges or centers. This was also corroborated by numer- 
ical simulations. After each DD cycle, we measured the 
remaining spin polarization by recording the NMR signal 
[HJ . During the signal acquisition, we applied continuous 
wave decoupling to prolong the FID and thereby increase 
the detection sensitivity. 

The decay of the signal amplitude as a function of the 
DD evolution time t = Mt c represents the survival prob- 
ability of the state p : 



s(t) = 



Tr{/W)(0} 
Tr{po/5 } 



where p(t) = U(t)poU'(i) is the density operator of the 
spin-system at time t. In the case of free evolution (no 
DD), this corresponds to the free induction decay signal 
represented in Fig. [5l The solid blue lines in panels (a) 
and (b) represent the FID signal of the 1 H and 13 C spins, 
respectively, which correspond to the observables I x and 
S x ■ 



(a) Longitudinal (b) Transverse 





Figure 4: (Color online) Different initial states of the spin 
system: The blue (thin, long) arrows indicate the orientation 
of the initial 13 C-spin polarization, and the red (thick, short) 
arrows show the direction in which the pulses are applied. 
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Figure 5: (Color online) Experimental FID of (a) H and (b) 
13 C spins without (blue solid line) and with (red dashed line) 
BLEW- 12 1 H homonuclear decoupling. The correlation time 
of the bath ( X H) increases by an order of magnitude if the 
BLEW-12 sequence is applied. The qubit system (the 13 C 
nuclei) remains almost unaffected by the BLEW-12 sequence. 
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B. Single parameter three pulse family 

We first consider the simplest sequence that contains 
unequal spacings between the pulses. Fig. [5] (a) and (b) 
show the time evolution of the 13 C magnetization from a 



Figure 6: (Color online) Time evolution of the C magneti- 
zation of the initial state (proportional to its survival proba- 
bility), for the (a) longitudinal and (b) transverse initial con- 
ditions subjected to UDD3. The value of -1 in the latter cor- 
responds to the state that is antiparallel to the initial state. 
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longitudinal initial state po = S x and a transverse initial 
state po = S y for different cycle times. They show that 
when the longitudinal case is compared to the FID, dy- 
namical decoupling allows for a longer survival of coher- 
ence, and this improves for shorter cycle times [HI. The 
signal decay is an order of magnitude faster for the trans- 
verse initial condition in Fig. [H^b), than the longitudinal 
condition [HI a) [38J. As we demonstrated in Ref. [38j], in 
this case the decay is due mostly to pulse imperfections 
(flip-angle errors), which lead to a loss of polarization. 
In the longitudinal case, the effect of these pulse errors is 
partly compensated over the sequence, but in the trans- 
verse case, they accumulate. This accounts for the large 
difference between the two cases already demonstrated in 
previous works [33L f38| . In particular in our experimen- 
tal setup, the flip-angle error was estimated to be around 

10% m. 

Uhrig introduced the idea that it might be advanta- 
geous to omit the condition that all pulse separations 
should be identical and derived a simple formula that 
determines the optimal spacing of N pulses. We exam- 
ine this scheme first in the simplest possible case, cor- 
responding to a cycle of 3 pulses. In order to determine 
the optimal pulse spacing experimentally, we parametrize 
the pulse spacings of the symmetric three-pulse sequence 
with a parameter shown in Fig. [TJa). x is defined as 
the deviation (as a fraction of the cycle time) of the pulse 
separation from those of the UDD3 sequence, which are 
Ti = T c sin 2 (7r/8) and T2 = t c /2 — t\. This parametriza- 
tion captures all possible mirror symmetric three-pulse 
sequences. When x = 0, the sequence is UDD 3 , and 
when x ~ 0, 0203, the sequence is CPMG. For this fam- 
ily of DD sequences, we experimentally determined the 
decay rate as a function of x by fitting a pure exponential 
function. 

Figure [8] shows the experimentally obtained decay rates 
for the longitudinal initial condition, as a function of 
the parameter x. We determined the value of x that 
generates the slowest decay rate by fitting the experi- 
mental data points with a quadratic function (not shown 
in the figure) and determining the position of the mini- 
mum. These minima are the red circles in Fig. [5] The 
sequence that generates the slowest decay is obtained 
for x m in sa 0.021 ± 0.002, for all cycle times used in 
the experiment. This is astonishingly close to the value 
xcpmg = 0.0203 that corresponds to the equally spaced 
sequence. 

In addition to the experimental data points, Fig. [5]also 
shows solid lines. They represent the results of numeri- 
cal simulations of the decay rates for the same parameters 
as the experiments. In the simulations, we assumed that 
the pulses were ideal. Following Ref. [U, [32j this in- 
volves modeling the pulse sequence as a one-dimensional 
filter /(t, t c ), as shown in[7f b). A detailed description of 
this model and the simulations is presented in Sec. I VII 
Since this model does not take into account pulse errors, 
it cannot be applied to the transverse initial condition, 
where these errors play a predominant role [H[. 
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Figure 7: (Color online) (a) Single parameter sequence family 
containing all mirror symmetric three-pulse sequences. They 
are described by the parameter x, where Ti and T2 are the ex- 
act UDD3 delays, (b) Time domain filter function fsir' ,2t c ) 
for two cycles of the UDD3 sequence, (c) Time domain filter 
function fi{r', r c ) for one cycle of the UDD4 sequence. In the 
filter representation we assume perfect pulses in the center of 
the finite-width pulses. The dashed lines indicate that the 
time domain filter toggles between ±1 at the position of the 
pulses. 
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Figure 8: (Color online) Experimental (symbols) and simu- 
lated (lines) decay rates of the magnetization as a function 
of the deviation x for different cycle times for the longi- 
tudinal initial condition. The UDD3 sequence corresponds 
to iudd = 0, while the CPMG sequence corresponds to 
icpmg ~ 0.0203. Note that the quadratic fitting curves used 
to determine to position of the minimum of each curve are 
not shown. 



For the transverse initial condition, the rates are about 
an order of magnitude larger than for the longitudinal 
initial condition, and they increase for short cycle times 
[38j . They are essentially independent of the parameter x 
because they are dominated by pulse errors. We therefore 
consider only longitudinal states in the following. 
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Determining the best three pulse DD sequence 



(a) Increasing power 
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Figure 9: (Color online) In order to determine the best longi- 
tudinal three-pulse DD sequence, the decay curves in Fig. [8] 
were fitted to a quadratic function in the region of the mini- 
mum. The positions of the minima of these quadratic fits are 
represented by the red circles in this figure. The black dotted 
line corresponds to icpmg ~ 0.0203. 



C. Variation of UDD order 

The performance of a dynamical decoupling sequence 
increases with the average power of the external control 
field used for DD Q . This means that given a fixed time 
window i, the DD efficiency improves with the number 
of refocusing pulses. Since the pulses are of finite length, 
the maximum number of pulses that can be employed 
during t is bounded [H M, El IH, EE HI • 

The scheme depicted in Fig. [TUTa) describes a situation 
in which the cycle time r c is fixed, while the UDD order 
is increased from top to bottom. Given the fixed window 
t = t c , the power used for DD during t increases from 
top to bottom. The corresponding survival probability is 
shown in Fig. [TOtc) for different UDD orders (from UDD2 
to UDD 14) and with r c = 600 /is (not including the pulse 
lengths). Clearly, if the cycle time is fixed, higher order 
UDD sequences maintain the qubit coherence for longer 
times. However, this improved performance is obtained 
at the price of an increase in the average power applied 
to the system. For many applications, this is not pos- 
sible, and a more meaningful comparison is obtained by 
keeping the average power level fixed. Fig. [TUJb), shows 
the corresponding pulse sequences and Fig. UOf d) the re- 
sulting decays. For these experiments, the cycle time was 
scaled with the UDD order N as r c — N 110.4 /is. Under 
these conditions, the variation of the decay rate with the 
UDD order is relatively weak and higher orders perform 
slightly worse than UDD 2 (CPMG) . 

We quantified the decay rates for UDD orders 2 to 30 
for t c = N 110.4 /xs. The results are shown in the right 
panel of Fig. [TT] The left panel schematizes how the DD 
power is kept constant across orders - the cycle time is 
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Figure 10: (Color online) Top panels: Pulse sequence scheme 
keeping fixed the time window t as the UDD order N is in- 
creased. To keep the DD power constant across UDD or- 
ders, the cycle time is scaled with TV (b). Bottom panels: 
Experimental 13 C signal decay for different UDD orders for 
(c) a fixed cycle time r c = 600/is, (d) a scaled cycle time 
t c = 110 AN fis. 



scaled with the UDD order. Clearly, for a large span of 
UDD orders, the CPMG sequence performs better than 
the UDD. The results mean two things: first, that in 
the regime of our experiment, higher UDD orders pro- 
vide no advantage over lower UDD orders; and second, 
that in consequence the UDD protocol itself performs 
worse than the CPMG sequence with the same number 
of pulses. Although these experiments were carried out 
with the average delay between pulses of 110.4/is, the 
same qualitative behavior is also observed for different 
times between pulses. However, the difference between 
the decay rates of CPMG and UDD can be reduced by 
shortening the cycle time as we show below. 



D. Variation of cycle time 

In the limit of infinitesimally short cycle time and in- 
finitely strong pulses - the bang-bang regime one can 
maintain quantum coherence for arbitrarily long times. 
However, in any real experiment, the peak power as well 
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Figure 11: (Color online) Left panel: Pulse scheme manifest- 
ing the cycle time scaling with the UDD order N. Right panel: 
Experimental (symbols) and simulated (line) decay rates of 
the 13 C magnetization for different UDD orders (blue circles) 
compared with decay rates achieved with the CPMG sequence 
(red rhombus). The same number of pulses during a time win- 
dow t axe applied by scaling the cycle time as r c = 110.4iV /xs. 



as the average power applied to the system are limited 
and should be minimized. Accordingly, the duration of 
the pulses and the cycle time cannot be reduced below 
some minimal values [l9|, HH, [42], |43|. In this section, 
we examine the DD performance of a given sequence as 
a function of the cycle time. The left panel of Fig. [12] 
schematically describes the variation of the cycle time for 
a fixed UDD order N, in this case for N = 5. In the right 
panel of Fig. I12[ the points are experimentally obtained 
decay rates while the lines are simulation results using 
the filter model for ideal pulses detailed in Sec. [Vlj For 
long cycle times, the performance of equidistant pulses 
(CPMG) is better than UDD however as r c is reduced its 
performance improves considerably, although the UDD 
decay rates are indistinguishable from the CPMG in this 
regime. Similar results were obtained for similar exper- 
iments with different UDD orders. Fig. [T^] shows that 
in the regime of our experimental setup, UDD does not 
perform better than the CPMG for any r c . 



E. Reducing the fluctuation rate 

The decay rate of the survival probability under dy- 
namical decoupling depends on the ratio of the cycle time 
r c to the correlation time of the bath tb [H, EH ■ The 
UDD sequence was designed to provide optimal decou- 
pling performance for baths with long correlation times. 
In our system, the bath correlation time depends on the 
strength of the dipolar coup lings dij between the 1 H nu- 



clear spins of the bath [37|. One can effectively rescale 
dij, and hence the correlation time of the bath, by ap- 
plying decoupling sequences to the bath spins [541 ] . This 
allows us to access regimes where the cycle time r c is sig- 



Figure 12: (Color online) Left panel: Pulse sequence scheme 
showing the variation of its cycle time for a fixed UDD order 
N — 5. Right panel: Experimental (points) and simulated 
(lines) decay rates of the 13 C magnetization for UDD5 and 
UDD20 and its comparison with equidistant pulse (CPMG) 
sequence having the same number of pulses for the different 
UDD cycle times. 



nificantly shorter than the correlation time of the bath.. 
For this we employ the well-known homonuclear decou- 
pling sequence BLEW- 12 67|. As shown by the dashed 
lines of Fig. [5] application of this pulse sequence in- 
creases the bath-correlation time by about one order of 
magnitude: it was around tb ~ 100/is without BLEW-12 
and tb ~ 1000/is with BLEW-12 for a cycle time t c b = 
84/xs. To determine the effects on the qubit system, we 
measured its free induction decay (FID). We found that it 
remains almost unaffected by the ^ homonuclear decou- 
pling (Fig. [5b). The scaling due to the BLEW sequence 
[671 ] . 0.475, appears to be indistinguishable from the ef- 
fect of self-decoupling, which reduces the 13 C linewidth 
in the absence of 1 H- homonuclear decoupling [68j]. The 
experimental scheme of Fig. [2] was used again, the only 
change being that the homonuclear decoupling sequence 
is applied to the bath in parallel with the application of 
the DD sequence to the system. 

Figure[I]|]shows the pulse sequence used during the DD 
period (left panel). For different UDD orders, the signal 
decays were measured with and without the BLEW-12 
sequence (dashed box in the left panel) on the 1 H chan- 
nel. The experimental results in the right panel of Fig. 
Q1J] show that the increase of the bath correlation time by 
roughly one order of magnitude leads to significant im- 
provements in the performance of all DD sequences. The 
performance of UDD3 and CPMG is indistinguishable, 
while that of UDDg is lower. 

Similar results are obtained for longer and shorter cy- 
cle times. This is evidenced in Fig. I14f b). which shows 
the survival probability after a single cycle of the UDD 
and equidistant pulse (CPMG) sequences for different cy- 
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Figure 13: (Color online) Left panel: DD pulse sequence 
scheme for UDD6 and CPMG applied in conjunction with 
the BLEW-12 sequence on the H spin-bath. The cycle time 
of the DD sequences is 6r, while for the BLEW-12 sequence 
the cycle time is t c b, where 6r/r c s is an integer. Right panel: 
13 C signal decay for the DD sequences described in the legend. 



cle times and DD orders. The pulse sequence for this 
experiment is shown in Fig. I14f a) for two cycle times 
t c = 336/iS and r c = 672/is of the UDD 12 sequence. The 
results show that for large UDD orders like TV = 12, the 
CPMG sequence performs better, even when the fluc- 
tuations of the environment are slow. For N = 3 the 
performance of both sequences is comparable within ex- 
perimental errors. 
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VI. PHYSICAL INTERPRETATION OF THE 
EXPERIMENTAL RESULTS 



A. Semiclassical approximation 



Figure 14: (Color online) (a) UDD12 sequence applied with 
homonuclear decoupling of the H spins. Sequences for two 
different cycle times r c = 336^is = 4r c s and r c = 672/is = 
8t c b, where r c s = 84/^s is the cycle time of the BLEW-12 
homonuclear decoupling sequence, (b) Experimental survival 
probability (fidelity) of the initial state pi = S x for a single 
cycle with TV pulses and varying the cycle time r c with and 
without homonuclear decoupling of the spin-bath. 



In this section, we present a semiclassical model to ex- 
plain and interpret the experimental results of section fVl 
Since the system is well in the high-temperature limit, 
it is possible to use a semiclassical description instead 
of the fully quantum mechanical treatment j2Z|. Most 
of the mathematical description of this subsection was 
developed within the DD context on different systems 
[24l [32I and they can be obtained from standard semi- 
classical treatments [37| . Here we connect and reinter- 
pret them to describe our spin system. Starting from the 
quantum mechanical description of scction|TTl the effect of 
the environmcnt-Hamiltonian % e on the evolution of the 
system may be discussed in an interaction representation 
with respect to the evolution of the isolated environment: 



the system-environment Hamiltonian then becomes 



-m E t 



UsEe 



S z e 



-iH E t 



JU E t 



b SE S z e- mEt E z e mEt 



(13) 
(14) 



where E, = 



^2 j j^li represents an effective spin bath 

operator and bsE = vSi ^Si * ne coupling strength. 
Since He does not commute with %se, the effec- 
tive system-environment interaction becomes time- 
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dependent: the system experiences a coupling to the en- 
vironment that fluctuates. For a semi-classical treat- 
ment we trace over the bath variables, and replace 
b SE e^ Et E z e l1lEt of Eq. © by the stochastic func- 
tion bsEE z {t) representing a classical random field with a 
Gaussian distribution with zero average, (E z (t)} = 0, and 
the autocorrelation function (E z {t)E z {t + r)) = g (r) 
[37j . The spectral density of the system- bath interaction 
is the Fourier transform of g(r), 



S(u) 



1 



<2ti J- 



dTg{r)e~ 



(15) 



It describes the relative weight of the different frequency 
components of the SE interaction. 

Using this effective field E z (t), we write the semiclas- 
sical SE interaction Hamiltonian as 



n SE (t) = b SE E z (t)S z 



(16) 



Clearly, T-LsE{t) commutes with itself at all times. This 
allows us to calculate the survival probability 



Tr 



s(t) 

by integrating 

rt 

dtiU SE {ti) 





'Po 



Tr{/5 /5 } 



b SE S z I dt 1 E z (t l )=4>{t)S z 
'0 



The survival probability then becomes 



»(*) 



Tr {popo} 



(17) 



(18) 



(19) 



If the spin ensemble is initially polarized along the z- 
direction, po = S z , S z is a constant of motion. However 
if po = S x ,y> its survival probability is 

Tl . L-i<Kt)s z g^ ye i<Kt)s z s Xt \ 

s x,y(t) = " f — -1 = COS </>(£). 



Tr 



Taking the average over the random fluctuations 



(20) 



(21) 



where we have used cos0(i) = (e'^W + e'^W) /2 and 

the property {e tX ) = e^ x ^^ x )/ 2 for a Gaussian ran- 
dom variable X. For a simple interpretation of the decay, 
we use the exponential's argument 



x(t) = \ m)) = lbsE 



t r-t 

dh / dt 2 E z (h)E z (t 2 ) 
Jo 



SE\ 



dh 



dt 2 g(h-t 2 )fo(t u t)f Q (t 2 ,t), 

(22) 



where f (t',t) = Q(t')Q(t ~ t') with the Hcavisidc func- 
tion O(i') = 1 V t' > and zero otherwise. Using Eq. 
(fT5|) . we rewrite this convolution integral as 



x(t) 



/ 2tF|6 



SE\ 



duS(uj) \F (u),t)\' 



(23) 



where F (uj,t) is the Fourier transform of fo{t',t) \2A 



F (u,t) 



/2tt J- 



dt'f (t',t)e 



-iujt 1 



/2tt 



, ^/ 2 sin(^/2) (24) 



( W /2) 



and we have used the fact that the random field E z is not 
correlated with the time-domain filter function /. The 
decay function x(t) is t nus equal to the product of the 
spectral density S(uj) of the system-environment coupling 
and the filter transfer function Fo(io,t). 

Here, we have treated the case where no control pulses 
are applied, which corresponds to the free induction de- 
cay. Since both S(uj) and l-F^o- 1 , i)| 2 reach their max- 
imum at ui = 0, it is the low frequency environmental 
noise that has the highest contribution to the decay rate. 



B. Analogy between FID and single slit diffraction 

Since the decay function %(i) arises from the interfer- 
ence between the random fields E z (ti), E z (t 2 ) at different 
times, it is helpful to draw an analogy with interference 
effects in optics. This analogy is best seen with the help 
of Huygcns' principle, which allows us to associate an 
elementary wave with every point in space; here, we ob- 
serve interference between elementary waves generated at 
different points in time and weighted by the filter func- 
tion F(oj, t). In the case of free precession (the FID), the 
time-domain filter function fo(t',t) is constant over the 
interval t' — [0,i\. This is exactly analogous to the case 
of diffraction from a slit that extends from to t, and 
as in the optical case, we obtain a diffraction pattern 
oc I sin(a;)/a;| 2 . 

According to Eq. (|2~i]) , the width of the diffraction pat- 
tern is oc j and its amplitude oc t. For very short times, 
i.e., narrow slits, the corresponding diffraction pattern is 
broader than the width of the spectral density S(uj) and 
the integral (|23|) grows ~ t 2 . When the time t exceeds 
the correlation time tq of the system-environment inter- 
action, the width of slit broadens and the width of the 
filter function Fo(u,t) becomes narrower than the spec- 
tral density pattern S(uj). In the long time limit ( > tb, 
the filter function narrows to a delta function at u) = 
and the decay function becomes x(i) oc \bsE\ 2 S(0)t, cor- 
responding to an exponential decay. This result is equiv- 
alent to the one obtained by Fermi's golden rule and is 
valid until a power law decay arises |69l [70| . 
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C. Effect of pulses : interference 

The effect of DD pulses is a modulation of the time- 
domain filter function f(t\ t) and therefore of the transfer 
or filter function F(u>,t). As shown by Eq (|2li)) . a slow- 
down of the decay is achieved by minimizing the overlap 
between S(oj) and \F(uj 7 t)\ 2 [H In the typi- 

cal case that the environmental spectral density peaks at 
small frequencies, this implies that \F(co,i)\ should be 
close to for small frequencies uj. 

Let us now study the effect of ideal DD pulses. They 
generate reversals of lisE^t), so that the resulting Hamil- 
tonian still commutes with itself at all times. As de- 
scribed in the previous sections, N pulses are applied 
during the interval r c at positions tj — {t\,t2, ■ ■ ■ , ijv}, 
with t = and tpf + i = r c . Under this condition, the 
time-domain filter function /n(t', t c ) in Eq. (j2"2")l be- 
comes 

N 

f N (t>, t c ) = ]T ®(*' - tj)e(t j+1 - t'), (25) 

j=o 

where /jv(£'j T c) now switches between ±1 at the position 
of every pulse [24| . This is depicted in Fig. [7(c) for a 
cycle of UDD4. We define r c / as the period of the filter 
function /jv(£',7")- For even DD orders, T c f is equal to 
the cycle time r c , for odd DD orders, it equals two cycles, 
i.e., T cf =2t c [Fig. mh)}. 

If the average of the function /iv(i', T cf) over the period 
T c f vanishes, the filter function vanishes at ui = 

2 



1^(0, T cf ) I 



1 



dt'f N (t',T c f] 



= 0. 



(26) 



For a general A^- puls e sequence, the transfer function 
F N (uj,t) becomes ]2~H] 



F N (uj,T c f) 



C/ / JV (t',r c/ ) e -^'dt' 



2tt Jo 



e~ luJt 'dt' 



1 i + (-i) Ar + i e-' : ^ + 2s;; 1 (-ii j 



j -iuitj 



2tt ioj 

(27) 

Uhrig found a suitable distribution of N pulses over 
a time T c f that eliminates the first N derivatives of 
u>F^{uj,T c f) [2^, HH and thus it has an optimally flat 
stop-band at uj = 0. 

For a given distribution of pulses, the corresponding 
time domain filter function is 

f N (t',T cf ) = Mi', 00) x Q{t')Q{T c f-t'), (28) 

where /jv(^', 00) is an infinite extension of the filter func- 
tion. /jv(t', 00) can be written as a Fourier series 

00 

f N (t',oc)= Y A k exp(iku> t'), (29) 

k— — 00 



where kuio = 27rfc/r c / are the harmonic frequencies of the 
period t c /, and 



Ah 



/2tt 

Tcf 



F N (ku ,T c f) 



(30) 



the amplitudes. Hence, by convolution, the frequency 
domain filter has the form, 

F N (uj,T c f) = 



Fn (fa*, ^f) e -^-^o)^/ 2 sii ;^- fc ^)^/ 2 ] 



(uj - ku ) T c f /2 



(31) 



To determine the behavior of F^(uj, r c /) close to u> = 0, 
it is sufficient to consider the effect of only the first few 
harmonics on cither side of w = 0. This is shown in 
Fig. rT5]for UDD4. The red thick line in the lower panel 
illustrates the stop-band for UDD4. The blue bars are 
the Fourier coefficients A-^ of the filter function f4(t\ 00). 
The contributions of the individual terms in Eq. pip arc 
represented by the dotted lines in Fig. [15] The sum 
(interference) of these diffraction effects gives rise to a 
maximally flat filter shape \Fn(uj 1 t c )\ close to w = 0. 
The top panel of Fig. rj5] shows the effect of the addition 
of diffraction effects of the first two harmonics to the right 
of uj = 0. The resulting (red thick line) in a region close 
to uj = is canceled by the corresponding harmonics to 
the left of uj = 0, as shown in the lower panel. 

Physically therefore, the sum in Eq. (f2"Tf can be under- 
stood as the interference of the diffraction effects due to 
each intcr-pulsc delay. In an optical analog, the switching 
between fN{t',T c f) = 1 — > f^(t',T c f) = —1 corresponds 
to a phase shift by tt, which could be implemented by 
a series of A/2 retardation plates. For short cycle time 
T c f <^ tb , the width of this region, where the filter func- 
tion vanishes, becomes broad compared to the S(uj) and 
the integral (|2"3"|) and thus the decay tend to zero. 



D. Effect of cycle iteration: diffraction grating 

Let us consider a DD sequence iterated M times for a 
total period tm = MT c f. This is illustrated in Fig. [T] 
The corresponding filter function is then 



F n (uj,t m = MT c f) 



1 



Mt',T c/ ) e - iW *'dt'x 



2tt Jo 

I + p -i^T c _|_ e -i2wr c/ _|_..._|_ e -i(M-l)uiT ef 



Mur cf 

sin ( — 2 — j iui{M~i)T cf \ 
-e 



sin I 



'2tt Jo 



/iv(t',T c/ )e- Mi '^' 



/ MU)T f \ 

3111 [—2^ ) 



sin I 



™(M-l)r cf 

e 5 F N {uj,T cf ). (32) 
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Figure 15: (Color online) Decomposition of the filter func- 
tion Fjv(w,t c /) for the UDD4 sequence. Blue bars mark the 
contributions of the harmonics of uo = 27r/r c /. Orange and 
magenta dotted lines labeled a and b are sine functions cen- 
tered at ujo and 2cjo- The top panel considers the effect only 
of the first two harmonics (k = 1,2). They interfere to give 
the resulting red thick line. The bottom panel shows how a 
region near to to = is completely canceled when we also con- 
sider the contribution of the two harmonics with k = — 1, —2. 
The green dashed line is a typical Gaussian spectral density 
function of the spin-bath. 



This is analogous to the intensity pattern obtained due 
to diffraction from an M-linc grating. The maxima of 
the filter function \Fn(uj,tm)\ again occur at the har- 
monic frequencies uj = 27rfc/r c /, where k is an integer. 
These are the peaks of the blue bars in Fig. [TS] and 
Fig. HSJa) for UDD 4 , and the red bars in Fig. \W[d) for 
CPMG. Their amplitudes are given by the filter function 
-F/v(w,t c /) of a single cycle which is shown by blue and 
red dotted lines in Fig. [ToTa) and (b). For different even 
UDD orders the single cycle filter functions Fjv(w,r c /) 
are shown in Fig. 116( b) and the respective harmonic po- 
sitions are represented by the empty circles. Between 
two of the principal maxima (harmonics) are (M — 2) 
secondary maxima, which are determined by the grating 
transfer function sin(Ma>T c //2)/sin(a;T c //2). As shown 
in Fig. [ToTa) and (d) on the solid lines and in panel 
(c), their amplitudes with respect to the principal max- 
imum falls off oc M _1 . For example, with only 6 cycles, 
the intensity of the first secondary maximum is less than 
5% of the maximum. Note that in Fig. [16] \Fjf((j,T c f)\ 
is plotted; the contribution of the secondary maxima is 
drastically reduced after taking the square. 

Hence for even a few cycles, the filter function 
\Fn(uj, tm)\ becomes an almost discrete spectrum that 
is given by the discrete Fourier transform of /iv(i')T c /). 
This is just the function /jv(i',oo). Fig. [TBI shows the 



comparison of the filter functions Fn(uj,tm) for UDD4 
[panel (a): blue solid line] with M = 12 and CPMG 
[panel (d): red solid line] with M = 24. The latter choice 
allows us to have the same evolution time for a fair com- 
parison. The luq in the frequency axis is defined in terms 
of the cycle time for CPMG, 



ujo = 2n/(r? f PMG ) = 2 x 2n/(r^) 



UDD4 



(33) 



Figure [ToTc) shows Fn(lj 7 tm) for different UDD orders 
iterated to match the same total evolution time. In Fig. 
[T5] and Fig. 116( a) and (d), the spectrum corresponding 
to F)v(w,oo) contains only the blue bars (for UDD4) or 
the red bars (for CPMG). Equivalently on Fig. IToTb). it 
is given by the empty circles. The decay function x(t) is 
thus a weighted sampling of the spectral density function 
S(oj), where the weighting factor is the magnitude of the 
respective Fourier components. Thus, 



X(t = t m ) oc t m b 2 3B ^2 A\S (kuj ) 



(34) 



Here, we have assumed t = tm S> tb so that F(ui,t) 
is well represented by a series of ^-functions centered at 
kujQ, given by Fn(lu, 00), while the contributions from the 
secondary maxima can be neglected. 

This decay function (fM]) is equivalent to the one de- 
rived by standard time dependent perturbation theory 
for a spin interacting with a continuum of states. The re- 
sult is similar to an expression derived by Fermi's golden 
rule but the spectral density S(ui) is evaluated on the fre- 
quency components of the time-dependent perturbation 



n SE (t) 



^ A k cos (ku t) 



E z {t)S z 



(35) 



It follows that for large M, the main contribution to the 
decay rate is determined by the spectral density S(lu) 
at the harmonic frequencies ktvo. The lowest frequency 
component present in a DD sequence is the inverse of the 
period r c f of the time domain filter. The CPMG sequence 
has the shortest period if the average distance between 
pulses r is fixed: r c ( } PMG = 2r and w CPMG = tt/(t). For 
an N pulse UDD cycle the period of the filter function in 
the time-domain is r c ^ DDjv = Nt (t c ^ DDjv = 2Nt) for N 
even (odd). Thus the first non- vanishing component in 
its Fourier expansion is at u^ BBn = 2tt/(Nt) (^ DDjv = 
it/{Nt)) [Fig. [TB]. In the system studied here, this has 
resulted in the CPMG generating the slowest decay. 



E. Simulations 

For the longitudinal initial condition flip- angle errors 
are well compensated as was shown in Ref. [38[. As a 
consequence, the effect of static pulse imperfections and 
finite pulse lengths in the experiments does not cause a 
qualitative difference to the simulations with ideal pulses. 
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Frequency w 



Figure 16: (Color online) Comparison of filter functions 
\F n (uj, t m )\ for different UDD orders: (a) UDD 4 with M = 1 
(dotted line) and M = 12 (solid line), even UDD orders -/V 
with M = 1 (b) and M = A&/N (c), and (d) CPMG=UDD 2 
with M — 1 (dotted line) and M = 24 (solid line), ojq is 
defined in terms of the cycle time of the CPMG sequence, 
ojq = 2tt/t^ mg . Blue circles (a), red rhombuses (d) and 
empty circles in (b) are the coefficients of the Fourier ex- 
pansion of /jv(£', oo). They are modulated by the shape of 
the filter function |_Fjv(ai, t c /)|, shown in panel (b) and rep- 
resented by blue and red dotted lines in panels (a) and (d), 
where rf/ MG = 2r and r c U ° DN = Nr. 



We therefore assume ideal DD pulses and Eqs. (|23j) and 
(f3"2j) allow us to simulate the resulting signal decays. 
However, for the transverse initial condition, the pulse 
errors are not compensated (38| . Thus, here we only de- 
scribe how the simulations for the longitudinal case were 
carried out. 

First, the spectral density 5(w) used in Eq. ([23")) is es- 
timated using the FIDs of Fig. [5J as described in section 
IV. A of Ref. [38| . We model a finite system to simulate 
the respective correlation functions that produce the ex- 
perimentally measured FIDs for both nuclei (see [38[ for 
details). The correlation time of the bath is defined as 



the time when the correlation function falls to 1/e of the 
initial value; this was found to be around tb = 110/xs. 

Now, including the DD pulses, x(i) is calculated after 
every cycle period r c following Eq. (|23")) and (|32p. The 
decay rate is determined by linearly fitting the values of 
x(t) obtained for different times. This procedure is exact 
for the case when the decay is purely exponential - for 
example, when the time t far exceeds the correlation time 
of the bath as described above. These simulations (lines) 
show reasonable agreement with the real experimental 
results in Figs. H HU and 



F. Discussion 

Visualizing dynamical decoupling as a filter [24|, HH 
for different frequency components of the environmen- 
tal noise provides a useful means for predicting relative 
DD performance. Some recent discussions by Biercuk 
and Uys 25| about the filter properties of different DD 
sequences concerning their Fourier components has ap- 
peared in parallel with our work. However, their work 
and ours contribute from different approaches: we com- 
pare the performance of different DD orders under the 
condition that the number of pulses applied during a time 
interval remains constant, while Biercuk and Uys fix the 
cycle time for every DD order. 

Although the UDD sequence that was derived from 
Eq. (|27|) has a flatter band-stop region close to ui = 0, 
it achieves this by the interference of a larger number 
of harmonic frequencies lower than the first harmonic 
of an equidistant sequence (CPMG) with the same av- 
erage spacing r between the pulses. This is illustrated 
in Fig. 1161 The first non-zero Fourier component of 



the CPMG occurs at w CPMG = tt/t l while for an TV- 
pulse UDD sequence it occurs at 



udd, 



,UDD od 



W CPMG , 



,CPMG 
' J 



/{2N)\. Therefore, 



/N 

as one increases 



the UDD order, additional components appear in the cor- 
responding frequency-domain filter function at frequen- 
cies below w CPMG . 

We may thus compare the effect of two contributions to 
the decay rate: The spectral density close to the lowest 
frequency component ujq and the integral over the fre- 
quency band close to ui = 0. UDD is designed to outper- 
form CPMG in this low-frequency band and thus superior 
if this is the dominant contribution. However, in the sys- 
tem that we arc considering here, the spectral density of 
the environmental noise is still sufficiently large at u1q DD 
that this term dominates and leads to UDD performing 
worse than CPMG. 

Previous works [ID, |43| used perturbation treatment 
of the SE interaction to predict the time evolution of 
the system under UDD at short times with high preci- 
sion. However, the regimes where the decay is reduced by 
increasing the UDD order lead in general to differences 
between their decays of order lower than 10 which are 
too small to be determined experimentally. The extra- 
polation of those decays to longer times, by means of 
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iterating the DD sequences, that may result in a mag- 
nification of the UDD performance relative to CPMG is 
not valid within the perturbativc treatment. Our anal- 
ysis, in contrast, is valid for times t 3> tb, where decay 
becomes appreciable (>1%) and thus experimentally ac- 
cessible. The filter function Fn(lu, Mr c f) now becomes 
discrete when M is large and changes drastically the de- 
cay compared to the extrapolation from a perturbativc 
treatment of the SE interaction. From a different ap- 
proach this situation was also observed by Khodjasteh 
et al. [45| for a bosonic bath: they noticed that an ex- 
trapolation could not be done and they obtained bounds 
for the UDD performance as a function of its order. Our 
results suggest that to achieve a parameter range where 
UDD would outperform CPMG would correspond to an 
almost static bath, i.e., the system-environment interac- 
tion |&se| 2 should be much stronger than the relevant 
intra-bath interactions. 

Our model may also explain the observation in (33j 
where the odd UDD orders were found to perform sub- 
stantially worse than even orders in a qubit system con- 
sisting of 9 Be + ions confined in a Penning trap, and where 
the S(u>) ~ 1/w 4 . This even-odd asymmetry was unex- 
plained in their paper [33j . 

VII. CONCLUSIONS 

We have presented experimental and theoretical results 
evaluating the relative performance of DD sequences with 
non-equidistant pulses in a purcly-dephasing spin-bath. 
We find that over a large range of cycle times, number of 
pulses per cycle and bath-correlation time, the equidis- 
tant sequence (CPMG) provides a measurably superior 
performance over non-equidistant sequences like UDD. If 
the rotation axis of the DD pulses is oriented in the di- 
rection of the initial state, we find that the experimental 
results match very well with simulations using the filter 
model for perfect pulses [24j|. In this case, the effects due 
to the finite pulse-length, and flip-angle errors do not 
qualitatively affect the results. 

We interpret our results through a semi-classical model 
by drawing an analogy of dynamical decoupling with in- 
terference effects in optics and filter theory. The CPMG 
sequence has the shortest cycle and its filter has the 



widest stop-band. By suitable interferences, the UDD 
filter has the flattest stop-band close to u = 0. In our 
system, the position of the lowest frequency component 
is more important than the very flat behavior of UDD, 
and thus the performance of CPMG is better than that of 
non-equidistant DD sequences. We show that this comes 
from the fact that the decay rate changes drastically from 
short to long times compared with the bath-correlation 
time. While for short evolution times, the flattest stop- 
band close to uj = may play an important role in re- 
ducing the decay rate, the fidelity reduction is < 1CT 4 
and thus too small to be observed experimentally. The 
extrapolation of those decays to longer times that may 
result in a magnification of the UDD performance rela- 
tive to CPMG is not valid. In contrast, we showed that 
for long times, where decay becomes appreciable, the fil- 
ter function becomes discrete, making CPMG superior 
because its frequency components have the widest spac- 
ing. Conversely, we expect that UDD should perform 
better in systems where the cycle time remains signifi- 
cantly shorter than the bath correlation time. 

Our results and the optical analogy we considered reaf- 
firm that it would be advantageous to tailor the DD 
sequence to the spectral density function of the noise 
in the particular system of interest. This is the basis 
behind the Locally Optimized DD (LODD) sequences 
II 11 [1>|, il III or similar [lEHi In particular LODD 
can be implemented using measurement feedback in or- 
der to find the optimal DD sequence (28|. However, more 
work needs to be done in methods for experimentally de- 
termining the exact noise features in the QIP system of 
choice, so that techniques like LODD could be used to 
determine optimal sequences for such environments. 
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